Dynamics of evaporative cooling in magnetically trapped atomic hydrogen 
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We study the evaporative cooling of magnetically trapped 
atomic hydrogen on the basis of the kinetic theory of a Bose 
gas. The dynamics of trapped atoms is described by the 
coupled differential equations, considering both the evapo- 
ration and dipolar spin relaxation processes. The numerical 
time-evolution calculations quantitatively agree with the re- 
cent experiment of Bose-Einstein condensation with atomic 
hydrogen. It is demonstrated that the balance between evap- 
orative cooling and heating due to dipolar relaxation limits 
the number of condensates to 9 x 10* and the correspond- 
ing condensate fraction to a small value of 4% as observed 
experimentally. 
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The recent development of atom-manipulation tech- 
niques [01 has realized Bose-Einstein condensation (BEC) 
in magnetically trapped alkali-metal atoms and 
atomic hydrogen In all these experiments, evapora- 
tive cooling was adopted at the final stage of the cooling 
procedures and was essential in obtaining the extremely 
low temperatures needed for the quantum degeneracy. 
The cooling mechanism of this powerful method is based 
on both the selective removal of energetic atoms through 
evaporation and coUisional rethermalizations among the 
remaining atoms Evaporative cooling itself is stim- 
ulating theoretical studies of the process of condensate 
formation in nonequilibrium atomic gases 

In BEC experiments with atomic hydrogen, evapora- 
tive cooling has been implemented by lowering the po- 
tential height of the saddle point at one end of a mag- 
netic trap with cylindrical symmetry This approach 
suffers from a reduction of the dimension of evapora- 
tion at low temperatures due to low elastic collision rate 
p8| , p^ since the s-wave scattering length, a, of a hydro- 
gen atom is anomalously small, about two orders of mag- 
nitude smaller than that in alkali- metal atoms. Heating 
caused by the dipolar spin relaxation then easily retards 
the evaporative cooling and prevents further cooling be- 
fore reaching the critical temperature of BEC. To over- 
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come this problem. Fried et al. adopted the evaporative 
cooling induced by radio-frequency (rf) magnetic field 
just after the "saddle-point" evaporative cooling The 
atom ejection technique utilizing the transition between 
trapped and untrapped spin states enables efficient three- 
dimensional evaporation even at low temperatures and in 
a highly anisotropic magnetic potential j^]. BEC was fi- 
nally achieved as a result of the heating-cooling balance 
in this efficient "rf-induced" evaporative cooling. 

There are several characteristic features of BEC in 
atomic hydrogen as compared with that in alkali-metal 
atoms 11. The small mass, m, of hydrogen resulted in the 
transition temperature of 50 /iK, the highest among BEC 
experiments. While the condensate fraction, i.e., the ra- 
tio of condensates to the whole trapped atoms, remained 
a small value of about 5%, a huge number of condensates 
containing about 10^ atoms were observed and the peak 
density of the condensates was about 25 times higher 
than that of the noncondensed atoms. Furthermore, with 
respect to the loss of condensates, both the three-body 
recombination process and the background gas collisions 
in a cryogenic environment of trap are negligible in the 
hydrogen system. Only dipolar spin relaxation becomes 
the dominant loss mechanism, and causes serious heating 

In this paper, the kinetic theory for evaporative cool- 
ing of a magnetically trapped Bose gas ||T^ is extended to 
include dipolar relaxation loss. This theory enables the 
consistent and quantitative investigations on the whole 
evaporative cooling process of the experiments, from the 
classical regime at the early stage of cooling to the quan- 
tum degenerate regime after the BEC transition. We ap- 
plied it to the rf-induced evaporative cooling in the recent 
BEC experiment with atomic hydrogen and quantita- 
tively investigated the characteristic features of hydro- 
gen BEC mentioned above. We expect further that our 
calculation studying the dynamics of trapped atoms in 
cooling process will give a good estimation for the grow- 
ing of condensates in the experiment, since there is room 
for discussions on the number of condensates measured 
in the spectroscopic way jS) . 

Here, we briefly mention the formulation of the theory 
. During evaporative cooling, the applied rf-magnetic 
field effectively truncates the trapping potential, U{r), 
at the energy, et, determined by the rf-field frequency. 
The thermalized distribution of noncondensed atoms in 
such a truncated potential is well approximated by the 
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truncated Bose-Einstein distribution function 
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e-i(r)-exp(ep//cBr)-l 



e(A(r)-ep), (1) 



where T is the temperature, Ep = /2m is the kinetic 
energy of atoms, and 6(x) is the step function, ^(r) 
denotes the local fugacity including the mean field po- 
tential energy; ^(r) = exp{[/^ — L/(r) — 2wn(r)]/fcsr}, 
where fx is the chemical potential, v = AiraTi^ /m is in- 
teraction strength between trapped atoms, and n{Y) is 
the density profile of atoms in a truncated potential. 
The step function in Eq. eliminates the momentum 
states whose kinetic energies exceed the effective poten- 
tial height A{y) = et — U{r) — 2vh{r). The momen- 
tum integration of this truncated distribution function 
gives the density profile in a self-consistent way such that 
n(r) = 4:Trh~^ J f{r,p) p^dp, and similarly the internal 
energy density such that e(r) = 47r/i~^ J €pf{r,p) p'^dp + 
vn^{r) + U{r)n{r). The total number of atoms, N, and 
the total internal energy, in a truncated potential 
are evaluated respectively by the spatial integrations of 
these density functions, N = J ri(r)dr and E — J e(r)dr. 
At low temperatures, the density profile in the con- 
densed region is described as the sum of condensates, 
no(r), and saturated noncondensed atoms, fin, such that 
n(r) = no{r) +nn- Condensates obey the Thomas-Fermi 
distribution no (r) = np — U{r)/v, where rip represents the 
peak density of condensates at the center of the poten- 
tial. The density n„ is evaluated through the truncated 
Bose-Einstein distribution function under the condition 
|(r) = l. 

The dynamics of trapped atoms during evaporative 
cooling is investigated on the basis of kinetic theory of 
a Bose gas Hydrogen atoms are removed from the 
trapping potential through both evaporation and dipo- 
lar spin relaxation processes The change rates 
(i.e., loss rates) of density functions, riioss and eioss, are 
evaluated respectively as the sum of the contributions of 
both processes, such that 

"-ioss(r) = -"ov(r) - G2 ■ K{r) ■ n^(r), 

eioss(r) = -eev(r) - G2 ■ K{r) ■ g(r) • n(r), (2) 

where ricv and ecv are the evaporation rates of density 
functions, G2 is the dipolar decay rate constant, and K 
is the correlation function which describes the second- 
order coherence of trapped atoms pO|-p^ . Both ricv and 
Bev were derived from a general collision integral of a Bose 
gas system in Rcf. and we adopt here the opposite 
sign of the notations of these rates defined in this refer- 
ence. The bosonic feature of evaporation process is there- 
fore included as the strong enhancement of these evap- 
oration rates due to the stimulated scattering of atoms. 
Regarding the dipolar spin relaxation process, we assume 
the correlation function of an ideal Bose gas such that 
K{r) = {nlir) + 4no(r)n„ + 2nl)/2ii^{r) [|o|-|2|. The 



spatial integration of Eq. (g) gives important thermody- 
namic quantities such as the loss rate of total number 
of trapped atoms, Moss, and that of the total internal 
energy, E^ioss: 



(3) 



Moss = / ?^loss(r) dr = -Mv - A^dip, 



E, 



eioss(r) dr = -E'^v - -Edip, 



where Ncv{Ecv) and iVdip(Edip) represent the evapora- 
tion rate and dipolar relaxation loss rate, respectively. 
These loss rates are complicated functions of tempera- 
ture T, chemical potential /i, and truncation energy et- 
Finally, the system obeys the coupled differential equa- 
tions 
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loss ; 
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(4) 



Time-evolution calculations were performed numerically 
on the assumption of quick rethermalization, which is 
appropriate for the slow evaporative cooling used in the 
usual BEC experiments Thus, the system always 

stays in a quasi-thermal equilibrium state and is de- 
scribed by the truncated Bosc-Einstcin distribution func- 
tion in Eq. (|l|). 

Precise knowledge of experimental conditions is also 
necessary for our quantitative calculations. The mag- 
netic trap adopted in the BEC experiment with atomic 
hydrogen ||] is modeled by the loffe-Pritchard potential 
U{p, z) = ^(ap)2 + (/3z2 + ey - 9 with the radial po- 
tential gradient a = 2.2 x 10"^^ J/ni, the axial curva- 
ture 2/3 — 6.8 X 10^2'' J/m , and the bias energy 9 = 
4.8 X 10^2* J (i.e., 9/kB = 35 ^K). At low energies, this 
potential is well approximated by the harmonic one with 
the radial oscillation frequencies u)p = aj \J m9 = 27r x 3.9 
kHz and the axial frequency uj^ = \/2pfm = 27r x 10 Hz. 
It should be noted here that, even after the BEC transi- 
tion, many noncondensed hydrogen atoms still distribute 
in the high-energy region, where such harmonic approx- 
imation fails. We treated the loffe-Pritchard potential 
exactly for the quantitative description of noncondensed 
atoms. Accordingly, all density functions were calculated 
in the cylindrical coordinate. 

As is well known, the efficiency of evaporative cool- 
ing strongly depends on the way the rf-field frequency is 
swept (i.e., forced evaporative coohng) ||]. The high effi- 
ciency is obtained by the slow evaporative cooling which 
continues for the time on the order of 10 seconds. In the 
BEC experiment with atomic hydrogen, the rf-field fre- 
quency was swept according to the functional form |23] 



Ht) = ( - 



where m — 23 MHz is the initial rf- 



frequency, Vf = 2 MHz is the final one, r = t/to denotes 
the normalized time with the specific time tg = 25 s, and 
g = 1.5 is the fixed parameter. The truncation energy et 
then changes according to the relation et(t) — hv{t) — 9. 
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Furthermore, we employed the following experimental 
parameters §: s-wave scattering length a = 6.48 x 10^^ 
nm p^ , th e dipolar decay rate constant G2 = 1.1 x 10~^^ 
cm'^/s |1^, initial temperature T = 120 fiK, initial peak 
density fi{0) = 5.0 x 10^'^ cm~'^, and the initial number of 
trapped atoms in a truncated potential N = 1.16 x 10^^. 
The time-evolution of evaporative cooling was calculated 
for 25 s, i.e., the normalized time r changes from to 1. 
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FIG. 1. Time evolution of (a) the total number of trapped 
atoms A'', and (b) temperature T. The arrow indicates the 
point at which the BEG transition occurs. The heating due 
to the dipolar spin relaxation occurs after the BEG transition. 

In Figs. 1(a) and 1(b), we show the time evolution of 
the total number of trapped atoms TV and that of temper- 
ature T, respectively. N decreases monotonically from 
1.16 X 10^^ to 8.7 X 10^ during a 25-s evaporative cool- 
ing. The BEC transition occurs at 15.4 s {t — 0.62) as 
indicated by the arrow in the figure. The decrease of 
temperature (from 120 to 22 /iK) is not so large in com- 
parison with that has been demonstrated in alkali-metal 
atoms . We can see that the cooling speed strongly 

slows down after the BEC transition in Fig. 1(b), indi- 
cating considerable heating occurs. Here we show the 
calculated results at the BEC transition point: temper- 
ature T = 56 fjK, peak density n(0) = 2.0 x lO" cm'^, 
number of trapped atoms N — 3.3 x 10^", and trunca- 
tion energy et/fce = 302 fiK. These parameters are quite 
consistent with the experiment which confirms the 
validity of our numerical calculations. 
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FIG. 2. Time evolution of (a) the total number of trapped 
atoms A*', and (b) temperature T. The arrow indicates the 
point at which the BEG transition occurs. The heating due 
to the dipolar spin relaxation occurs after the BEG transition. 



The condensate formation process shows the character- 
istic feature depicted in Fig. 2. Condensates first grow 
very rapidly after the transition. We can explain this by 
the fact that bosonic stimulation in evaporative cooling 
process [|l^ strongly accelerates the growing speed, since 
the density of condensates becomes very high in atomic 
hydrogen. The number of condensates A^o reaches the 
maximum value at 17.2 s (r = 0.69), and then decreases 
gradually due to dipolar decay. The calculated parame- 
ters at this maximum point are: temperature T — A9 /iK, 
peak density of condensates rip = 4.3 x 10^^ cm^"^, peak 
density of noncondensed atoms fi„ ~ 1.7 x lO^'' cm^"^, to- 
tal number of condensates Nq = 9.0 x 10^, total number 
of noncondensed atoms iV„ = 2.0 x 10^°, and truncation 
energy et//cB = 239 fiK. The corresponding condensate 
fraction at the maximum point of Nq is calculated as 
/ = No/{No + Nn) to be a small value of 4.3%. These 
results quantitatively agree with the recent observation 
of BEC with atomic hydrogen |^ . 

The number of trapped atoms has been evaluated ex- 
perimentally by measuring the density-dependent fre- 
quency shift of the two-photon 15-25' transition in atomic 
hydrogen |l^. The condensate number, iVo ~ 10^ in Ref. 
Q, was obtained on the assumption that the frequency 
shift for condensates is as large as that for noncondensed 
atoms with the same density. If we expect the disappear- 
ance of exchange effects in the excitation of condensates 
as discussed in Ref. the frequency shift for conden- 
sates becomes half of that for noncondensed atoms, and 
many more condensates, A^o — 6 x 10^, with an unrea- 
sonably higher condensate fraction, / — 25%, are given. 
The calculated results in Fig. 2 support the evaluation of 
condensate number which assumes the equal frequency 
shift for condensates and noncondensed atoms. 

Next, we discuss the influence of dipolar loss on evapo- 
rative cooling. From Eq. (H), the dipolar relaxation pro- 
cess preferably removes atoms in the region with higher 
density. Since the removed atoms have lower internal 
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FIG. 3. Time evolution of (a) the dipolar loss rate A'dip, 
and (b) normalized rates, N^v/N (dashed line) and N^ip/N 
(solid line). The dipolar loss rate exceeds the evaporation rate 
at around 17.5 s. The arrow indicates the point at which the 
BEG transition occurs. 



energy in the trapping potential, the dipolar relaxation 
loss results in system heating. The time evolution of 
the dipolar loss rate A'dip is shown in Fig. 3(a). The 
curve exhibits the jump structure around the BEC tran- 
sition point indicated by the arrow. The serious heating 
after the transition in Fig 1(b) is caused by such dras- 
tic increase of dipolar loss known as "relaxation explo- 
sion" in hydrogen system pO| . To clearly demonstrate 
the competition between evaporation and dipolar relax- 
ation processes during evaporative cooling, the normal- 
ized rates Ng^/N and N^ip/N B-re plotted together in Fig. 
3(b). The evaporation rate is larger than the dipolar loss 
rate until just after the BEC transition, which indicates 
that evaporative cooling is sufficiently efficient in this re- 
gion. The explosively enhanced dipolar-loss rate after the 
transition finally exceeds the evaporation rate at around 
17.5 s, and the resultant heating becomes serious. We 
note here that this time almost corresponds to the point 
where the number of condensates shows the maximum 
value in Fig. 2. The heating-cooling balance therefore 
limits the producible number of condensates and the cor- 
responding condensate fraction. One can expect that the 
trapping potential with weaker confinement and the op- 
timized evaporative cooling will move this balance in the 
direction of a larger condensate production p3|] . 

In conclusion, we have investigated the rf-induced 



evaporative cooling of magnetically trapped atomic hy- 
drogen on the basis of the kinetic theory of a Bose gas. 
This approach can be applied over the whole temperature 
region of evaporative cooling while the earlier investiga- 
tions based on the classical kinetic theory |^-||,|l^,^ fail 
around the critical point of BEC transition. The calcu- 
lated results quantitatively agree with the recent BEC 
experiment of atomic hydrogen. This proves our calcula- 
tions very useful, and the important future work would be 
to optimize the cooling trajectory for a larger production 
of condensates Our theory, on the other hand, as- 
sumes that the system always stays in the quasi-thermal 
equilibrium states during evaporative cooling. A precise 
study of the deviation from the quasi-thermal equilib- 
rium would give the limit of the approximation in 
the theory. 
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